February 4, 2008 10:26 



WSPC - Proceedings Trim Size: 9in x 6in hoerande! 



1 



COSMIC RAYS FROM THE KNEE TO THE SECOND 
KNEE: 10" TO 10^^ eV * 

JORG R. HORANDEL 

University of Karlsruhe, Institute for Experimental Nuclear Physics, P.O. 3640 
76021 Karlsruhe, Germany — www-ik.fzk.de/r^joerg 



The energies of cosmic rays, fully ionized charged nuclei, extend over a wide 
range up to 10^" eV. A particularly interesting energy region spans from 10^* 
to 10^* eV, where the all-particle energy spectrum exhibits two interesting 
structures, the 'knee' and the 'second knee'. An explanation of these features 
is thought to be an important step in understanding of the origin of the high- 
energy particles. Recent results of air shower experiments in this region are 
discussed. Special attention is drawn to explain the principle of air shower 
measurements — a simple Heitler model of (hadronic) air showers is developed. 
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1. Introduction 

The energy spectrum of cosmic rays (fully ionized atomic nuclei) spans 
a wide range in energy from GeV energies up to lO^*' eV. Over these 10 
decades the flux decreases by about 30 orders of magnitude rather feature- 
less, following roughly a power law dN/dE cx E"' . The power law behavior 
indicates a non-thermal origin of the particles. To reveal small structures 
in the shape of the energy spectrum the flux is usually multiplied by the 
energy to some power. The energy spectrum multiplied by E^ is depicted in 
Fig. 1. In this representation the spectrum looks rather flat and fine struc- 
tures can be recognized, indicating small changes in the spectral index 7. 
The most important are the knee at Ek ~ 4.5 PeV where the power law 
spectral index changes from 7 = —2.7 at low energies to 7 « —3.1, the 
second knee at E2nd ~ 400 PeV« 92 x Ek, where the spectrum exhibits 
a second steepening to 7 « —3.3, and the ankle at about 4 EeV, above 
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Fig. 1. All-particle energy spectrum of cosmic rays, the flux is uniltiijlicd by E'^, for 
references see [1]. The lines represent spectra for elemental groups (with nuclear charge 
numbers Z as indicated) according to the poly-gonato model [2] . The sum of all elements 
(galactic) and a presumably extragalactic component are shown as well. 



this energy the spectrum seems to flatten again to 7 » —2.7. To under- 
stand the origin of these structures is expected to be a key element in the 
understanding of the origin of cosmic rays (CRs). 

The lecture starts with a short overview on the physics of galactic cosmic 
rays (Sect. 2). Measurements in the energy region of interest are performed 
with air shower experiments, their principles arc outlined in Sect. 3 with a 
simple Heitler model. Finally, recent results of air shower experiments are 
reviewed in Sect. 4. 

2. Galactic Cosmic Rays and the Knee 
2.1. Sources 

At energies around 1 GeV/n all elements known from the periodic table 
with nuclear charge number Z from 1 to 92 have been found in CRs [2 4]. 
Overall, the abundance of elements in CRs is very similar to the abundance 
found in the solar system, which indicates that CRs are "regular matter" 
but accelerated to very high energies. This is emphasized by measurements 
of the CRIS experiment [5] which show that the abundances of particular 
isotopes in cosmic rays and in the solar system differ by less than 20%. 

The bulk of CRs is assumed to be accelerated in blast waves of supernova 
remnants (SNRs). This goes back to an idea by Baade and Zwicky who 
proposed SNRs as cosmic-rays sources due to energy balance considerations 
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[6]. They realized that the power necessary to sustain the cosmic-ray flux 
could be provided when a small fraction ~ 10% of the kinetical energy 
released in supernova explosions is converted into CRs. Fermi proposed a 
mechanism to accelerate particles with moving magnetic clouds [7] . This led 
to todays picture that the particles are accelerated at strong shock fronts 
in SNRs through first-order Fermi acceleration [8-12]. This theory predicts 
spectra at the sources following a power law dN/dE oc E~'^-^. 

Diffusive, first-order shock acceleration works by virtue of the fact that 
particles gain an amount of energy AE' oc at each cycle, when a cycle 
consists of a particle passing from the upstream (unshocked) region to the 
downstream region and back. At each cycle, there is a probability that the 
particle is lost downstream and docs not return to the shock. Higher energy 
particles are those that remain longer in the vicinity of the shock and have 
enough time to achieve the high energy. After a time T the maximum energy 
attained is Emax ~ ZepsBTVg, where /3s = Vg/c is the velocity of the shock. 
This results in an upper limit, assuming a minimal diffusion length equal 
to the Larmor radius of a particle of charge Ze in the magnetic fields B 
behind and ahead of the shock. Using typical values of Type II supernovae 
exploding in an average interstellar medium yields Emax ~ .^-100 TeV [13]. 
More recent estimates give a maximum energy up to one order of magnitude 
larger for some types of supernovae Emax Z ■ 5 PeV [14-16]. As the 
maximum energy depends on the charge Z, heavier nuclei (with larger Z) 
can be accelerated to higher energies. This leads to consecutive cut-offs of 
the energy spectra for individual elements proportional to their charge Z, 
starting with the proton component. 

This theory is strongly supported by recent measurements of the HESS 
experiment [17,18], observing TeV 7-rays from the shell type SNR RX 
J1713. 7-3946. For the first time, a SNR could be spatially resolved in 7- 
rays and spectra have been derived directly at a potential cosmic-ray source. 
The measurements yield a spectral index 7 = —2.19 ± 0.09 ± 0.15 for the 
observed 7-ray flux. The results are compatible with a nonlinear kinetic 
theory of cosmic-ray acceleration in supernova remnants and imply that 
this supernova remnant is an effective source of nuclear CRs, where about 
10% of the mechanical explosion energy are converted into nuclear CRs [19]. 

2.2. Propagation 

After acceleration, the particles propagate in a diffusive process through 

the Galaxy, being deflected many times by the randomly oriented magnetic 
fields (B f» 3 /uG). The nuclei are not confined to the galactic disc, they 
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propagate in the galactic halo as well. The diffuse 7-ray background, ex- 
tending well above the disc, detected by the EGRET experiment, exhibits 
a structure in the GeV region, which is interpreted as indication for the 
interaction of propagating CRs with interstellar matter [20]. The 7-rays 
are produced in inelastic hadronic reactions of CRs with the interstellar 
medium (ISM) via neutral pion decay p + ISM ^ n'^ ^ 77. The height 
of the propagation region in the halo has been estimated measuring the 
-•^^Bc/^Bc-ratio with the ISOMAX experiment to be a few kpc [21]. Deter- 
mining the abundance of radioactive nuclei, which decay on the way from 
the source to the Earth, allows to determine the residence time of CRs in 
the Galaxy. Measurements with the CRIS instrument yield about 15 • 10^ a 
for particles with GcV energies [22]. 

Information on the propagation pathlength of CRs is often derived from 
the measurement of the ratio of primary to secondary nuclei. The latter 
ar(^ produc;cd through spallation during propagation in the Galaxy. The 
energy dependence of the measured ratio is frequently explained in Leaky 
Box models by a decreasing pathlength of CRs in the Galaxy A(i?) = 
Ao{R/Ro)~^, with typical values Aq » 10 - 15 g/cm^, ^ « 0.5 - 0.6, and 
the rigidity i?o « 4 GV [23]. In this picture the spectra observed at Earth 
should be steeper as compared to the source, i.e. the spectral index 7 should 
be smaller by the value of 5. 

Energy spectra of individual elements have been measured up to ener- 
gies of about 10^'' eV by experiments above the atmosphere, the results 
being well compatible with power laws [2,4]. Due to spallation during the 
propagation process, the spectra of heavy elements are slightly flatter as 
compared to light nuclei [2,24], e.g. comparing protons 7^ = —2.71 ± 0.02 
to iron nuclei -jFe = —2.59 ± 0.06. 

The regular component of the galactic magnetic field will cause particles 
with charge Z to describe helical trajectories with a Larmor radius Rl = 
p/{ZeBo) — 1.08 pc-i?[PeV]/(Z-i3o[AiG]), while the random field component 
causes diffusive propagation. With increasing energy (or momentum) it 
becomes more difficult to magnetically confine the particles to the Galaxy. 
Since Rl oc 1/Z it is expected that leakage from the Galaxy occurs for light 
elements (low Z) earlier as compared to heavy nuclei, i.e. protons leak first 
and subsequently all other elements start to escape from the Galaxy. 

2.3. Structures in the Energy Spectrum 

Many possible origins for the knee arc discussed in the literature [25,26]. 
Most popular are assumptions of a finite energy attained during the accel- 
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oration process and leakage from the Galaxy as discussed. In both scenarios 
the energy spectra of elements exhibit a cut-off at an energy proportional 
to the nuclear charge Z and the knee in the all-particle spectrum is caused 
by the cut-off of protons. 

All other elements follow subsequently and above a certain energy no 
more particles are left. On the other hand, the measured all-particle flux 
extends up to 10^° eV, and the highest-energy particles are usually being 
considered of cxtragalactic origin. The Larmor radius of a proton with an 
energy of 10^° eV in the galactic magnetic field is i?L « 36 kpc, compara- 
ble to the diameter of the Galaxy. This emphasizes that such high-energy 
particles arc of cxtragalactic origin. The transition region from galactic to 
cxtragalactic CRs is of particular interest, key features are the origin of the 
second knee and the ankle. 

Reviewing the properties of GRs accelerated in SNRs, Hillas finds that 
a second (galactic) component is necessary to explain the observed flux at 
energies above 10^^ eV [27]. Another possibility is a significant contribu- 
tion of ultra-heavy elements (heavier than iron) to the all-particle flux at 
energies around 400 PeV [2,24]. In this approach the second knee is caused 
by the fall-off of the heaviest elements with Z up to 92. It is remarkable 
that the second knee occurs at Eind ~ 92 x E'fc, the latter being the energy 
of the first knee. The dip seen in the spectrum between 10^^ and 10^^ eV, 
see Fig. 1, is proposed to be caused by electron-positron pair production of 
CRs on cosmic microwave background photons [28] p -|- 73^ — » p -t- e+ -|- e~ . 

3. Measurement Techniques 

To clarify the situation and to distinguish between the different models, 
measurements of the flux of individual elements, or at least groups of ele- 
ments, up to high energies are necessary. Direct measurements above the 
atmosphere on stratospheric balloons up to energies exceeding 10^^ eV 
are performed with various instruments Uke ATIC [29], CREAM [30], 
BESS [31], or TRACER [32]. The presently largest experiment with single- 
element resolution, TRACER, has an aperture of 5 m^ sr. With an exposure 
of 50 m^ sr d accumulated during a circumpolar flight in 2003 energy spec- 
tra could be measured up to about 5 • 10^^ eV for oxygen and 8 • 10^^ eV 
for iron nuclei [33]. 

To extend the measurements to energies beyond the knee, at present, 
ground based installations are the only possibility. With these experiments, 
secondary products generated in the atmosphere are measured, the ex- 
tensive air showers (EAS). Air showers were discovered in 1938 by W. 
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Kolhorster [34] and independently by P. Auger [35]. Auger describes his 
work in a book [36] translated by the director of this school M.M. Shapiro. 

In todays experiments, the energy is basically derived from the number 
of particles observed and the primary's mass is estimated by measurements 
of the depth of the shower maximum (Sect. 3.1.5) or the electron-to-muon- 
ratio (Sect. 3.1.6). Two types of experiments may be distinguished: installa- 
tions measuring the longitudinal development of showers in the atmosphere 
and apparatus measuring the density (and energy) of secondary particles 
at ground level. 

An example for the latter is the KASCADE experiment [37], covering 
an area of 200 x 200 m^ . The basic idea is to measure the electromagnetic 
component in an array of unshielded scintillation detectors and the muons 
in scintillation counters shielded by a lead and iron absorber, while the 
hadronic component is measured in a large calorimeter [38] . The total num- 
ber of particles at observation level is obtained through the measurement of 
particle densities and the integration of the lateral density distribution [39] . 
The direction of air showers is reconstructed through the measurement of 
the arrival time of the shower particles in the individual detectors. 

The depth of the shower maximum is measured in two ways. Light- 
integrating Cerenkov detectors hke the BLANCA [40] or TUNKA [41] ex- 
periments are in principle arrays of photomultiplier tubes with light collec- 
tion cones looking upwards in the night sky, measuring the lateral distribu- 
tion of Cerenkov light at ground level. The depth of the shower maximum 
and the shower energy is derived from these observations. Imaging tele- 
scopes as in the HiRes [42] or AUGER [43] experiments observe an image 
of the shower on the sky through measurement of fluorescence light, emitted 
by nitrogen molecules, which had been excited by air shower particles. 

3.1. A Heitler Model for Air Showers 

The basic properties of EAS are illustrated using a Heitler model [44], 
expanding an approach by Matthews [45] . The principle ideas of the model 
are emphasized by full EAS simulations using the CORSIKA code [46] with 
the hadronic interaction models FLUKA [47] and QGSJET 01 [48]. For the 
latter, a modification with lower cross-sections has been used [49]. 



''Vertical showers with fixed energies between 10''' and 3.16 ■ 10^" GeV in steps of half 
a decade have been calculated. Thresholds for photons, electrons, muons, and hadrons 
were chosen as E-y > 0.25 MeV, Ee > 0.25 MeV, > 100 MeV, and E^ > 100 MeV. 
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Fig. 2. Left: Schematic view of an electromagnetic cascade {left) and a hadronic shower 
(right). Not all pion lines are shown [45]. Right: Number of electrons at shower maximum 
and depth of the shower maximum as function of photon energy. The lines are according 
to (1) and (2). 

3.1.1. Electromagnetic Cascades 

A simple approximation of an electromagnetic cascade is shown schemati- 
cally in Fig. 2. A primary photon generates an e+e~ pair. An electron radi- 
ates a single photon after traveling one splitting length d = Xq In 2, where 
Xq is the radiation length (Xg"" — 36.66 g/cm^). An electron looses on av- 
erage half of its energy through radiation over the distance d. After traveling 
the same distance a photon splits into an e'^e~' pair. In either instance, the 
energy of a particle is assumed to be equally divided between two outgoing 
particles. After n splitting lengths, at a distance x = nXo\n2, the total 
shower size (electrons and photons) is iV = 2" = exp(a;/A'o) and the initial 
energy Eq is distributed over TV particles. The splitting continues until the 
energy per particle Eq /N is too low for pair production or bremsstrahlung. 
Heitler takes this energy to be the critical energy {E^ — 85 MeV in air), at 
which ionization losses and radiative losses are equal. 

A shower initiated by a primary photon reaches its maximum size N„iax 
when all particles have the energy E'^, which means Eq = E^Nmax- The 
penetration depth X^ax at which the shower reaches its maximum is ob- 
tained by determining the number Uc of splitting lengths, required to reduce 
the energy per particle to E^. Since N^ax = 2"'', the number of splitting 
lengths is = \n{EQ/E'^)/\ii2, giving N,nax = E^/E^ and 

XI,, = n,X„ \n2 = Xo HEo/E^). (1) 

The elongation rate A specifies the increase of X^ax with energy Eq and is 
defined as A = dX,nax/d\gEo. Using (1) gives A^ = In 10 Xa = 84.4 g/cm^ 
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per decade of primary energy for electromagnetic showers in air. Thus, 
Xmax = 597 g/cm^ + 84 g/cm^ lg(£;o/PeV) is expected. This prediction 
agrees well with results of full simulations as can be inferred from Fig. 2. ^ 

The simple model describes quite well the position of the maximum of 
electromagnetic cascades when compared to EAS simulations and to mea- 
surements at accelerators [45,52]. However, the model overestimates the 
actual ratio of electrons to photons. It predicts that after a few genera- 
tions the electron size approaches ~ ^N^ax- This is much too large 
for several reasons, mainly that multiple photons are often radiated during 
bremsstrahlung and many electrons and positrons range out in the air. 

To extract the number of electrons iV™"^ at shower maximum from 
Heitler's total size N^ax, a simple correction 

is adopted, with a constant value g. When the estimated electron number is 
compared to measurements, the factor g has to be fine tuned. It depends on 
properties of the detectors used like the energy threshold and the efficiency 
to detect photons and electrons (or positrons) . Comparisons with results at 
accelerators indicate values between g = 10 [45] and g = 20 [52]. Results of 
a full EAS simulation are depicted in Fig. 2. For electromagnetic showers 
the number of electrons turns out to be almost exactly linearly proportional 
to the shower energy as expected from (2). A fit yields ex Eq-^'^ and a 
correction factor w 13 is obtained, compatible with the accelerator based 
results. With this value the number of electrons at shower maximum is 
according to (2) iV™"^^ = 9.0 • 10^ • Eo/PeY. 



3.1.2. Hadronic Showers 

Hadron induced showers are modeled using a similar approach, for a figura- 
tive sketch, see Fig. 2. The atmosphere is divided in layers of fixed thickness 
Ailn2, where Aj is the interaction length of strongly interacting particles. 
An energy around 100 GeV is a typical energy for pions in air showers and 
for a simple approach a constant value A, = 120 g/cm^ is adopted. Hadrons 
interact after traversing one layer, producing Nch charged pions and ^A^cfc 
neutral pions. The latter decay promptly to photons, initiating electromag- 
netic cascades. Charged pions continue through another layer and interact. 



''The deviations at high energies are due to the Landau Pomeranchuk Migdal effect 
[50,51]. 
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The process continues until the charged pious fall bellow the critical energy 
E^, where they are all assumed to decay, yielding muons. 

The multiplicity of charged particles produced in hadron interactions 
increases very slowly with laboratory energy oc in pp and pp data [53]. 
The multiplicity in tt-^^N collisions increases as Nch ~ 5, 11, and 27 at 10, 
100, and 10^ GeV, respectively [49]. A constant value Nch — 10 is adopted 
in the following for the number of charged particles produced in pion-air 
interactions, again corresponding to an energy of about 100 GeV. 

The second parameter is the energy at which further particle produc- 
tion by TT^ ceases. EJ! may be defined as the energy at which the probability 
for decay and hadronic interaction equalize. Following Rcf. [45] a constant 
critical pion energy EJ! = 20 GeV is adopted in the following. 

If we consider a proton with Eq entering the atmosphere, we have after n 
interactions N,^ = (Nch)"' charged pions. Assuming equal division of energy 
during particle production, these pions carry a total energy of (2/3)"£'o. 
The remainder of the energy goes into electromagnetic showers from tt*^ 
decays. Hence, the energy per charged pion is E^^ = Eo/{^Nch)"- After a 
certain number ric of generations, becomes less than E^. The number 
of interactions needed to reach E„ = E^ is 

ne = i^^ = 0.851gr^V (3) 
3.1.3. Number of Muons 

The number of muons is obtained, assuming that all pions decay, using 
A^^ = Nt^ = {Nch)""- Their energy dependence is derived applying (3) 

In AT^ = rie In Ar,^ = /3 In f^V with /? = « 0.85 (4) 

for Nch = 10. It should be noted that although Nch changes (slowly) as the 
shower develops, /? depends only logarithmically on this value. 

So far, an important aspect of hadronic interactions has been neglected. 
In an interaction only a fraction of the energy is available for secondary 
particle production, usually characterized by the the inelasticity k. Taking 
this effect into account, in an interaction initiated by a particle with energy 
E, the energy (1 — k)E is taken away by a single leading particle, ^kE is 
used to produce Nch charged pions, and ^kE goes via neutral pions into 
the electromagnetic component. Including inelasticity in the Heitler model 
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changes the parameter (3 in (4) to [45] 
ln[l + N,h] 



1 - 



K 



= 1 - 0.14k;. 



(5) 



In [(l + iV,;,)/(l-i«)] 



The elasticity for the most energetic meson in pion-air interactions yields 
1-K between 0.26 and 0.32 [49], resulting in /3 = 0.90. 

To expand the simple approach from primary protons to nuclei, the su- 
perposition model is used. A nucleus with atomic mass number A and en- 
ergy Eq is taken to be A individual single nucleons, each with energy Eq/A, 
and each acting independently. The resulting EAS is treated as the sum of 
A individual proton induced showers, all starting at the same point. The 
observable shower features are obtained by substituting the lower primary 
energy into the expressions derived for proton showers and summing A such 
showers. Applying this to the number of muons yields = A{Eo/{AE1!))^ . 
The number of muons in showers induced by nuclei with mass number A 
and energy Eq is then 



Two important features follow from (6): the number of muons increases 

as function of energy slightly less than exactly linear and iV^ increases 
as function of the mass of the primary particle as oc A"'^. Accordingly, 
iron induced showers contain about 1.5 times as many muons as proton 
showers with the same energy. This results from the less-than-linear growth 
of the number of muons with energy - (3 < 1 in (6). The lower energy 
nucleons which initiate the shower generate fewer interaction generations, 
and consequently, loose less energy to the electromagnetic component. 

The number of muons at shower maximum as function of energy is shown 
in Fig. 3 as obtained from full simulations. The lines indicate predictions 
according to (6), being well in agreement with the simulations. 

3.1.4. Number of Electrons 

Conservation of energy implies that the primary energy is split into elec- 
tromagnetic and hadronic parts Eq = i?em + E^- The number of electrons 
is estimated using this relation. The hadronic energy appears in the simple 
approach in the muon component as Eh = N^E^ and the energy fraction 
for the electromagnetic component is, using (6) 






^0 - N^El 

Eq 



= 1 - 
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Fig. 3. Number of muons (left) and number of electrons (right) at shower maximum 
as function of energy for primary protons and iron iniclei according to CORSIKA sim- 
ulations (symbols). The lines are predictions according to (6) and (10), respectively, for 
protons ( — ) and iron ( ) nuclei. 



The electromagnetic fraction is 57% at Eq = 10^^ cV, increasing to 79% at 
10^^ eV for proton induced showers. For iron induced showers the fraction 
rises from 38% to 68%. 

Equation (7) can be approximated by a power law 



Ee.m ( El 



Eo \AE^ 



(8) 



Series expansion near xq = Eq/E^ = 10^ yields the number of electrons at 
shower maximum as function of energy 

gE- gE-^ ^' ° ^' 

with6 = (l-/3)/(xo"''-l) « 0.046anda= (1 -a;^"^)/(a;(5) f« 0.40. Hence, 
a = 1 + b 1.046 is obtained, which leads, using g = 13 to 

Ne^5M-10'-A-°°^^(^^\ . (10) 

This implies that the number of electrons grows as function of energy 
slightly faster than exactly linear. The electron number decreases with in- 
creasing mass number, an iron induced shower is expected to contain about 
83% of the electromagnetic energy of a proton shower with the same energy. 
It should be emphasized that the model does not take into account absorp- 
tion in the atmosphere, thus, the number of electrons obtained is valid at 
shower maximum. 
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The number of electrons at shower maximum according to full simula- 
tions is shown as function of energy in Fig. 3. The results are compared to 
predictions according to (10) for proton and iron induced showers indicated 
by the lines. It can be seen that the simple model reproduces quite well the 
results of the full simulations. 



3.1.5. Depth of the Shower Maximum 

The atmospheric depth at which the electromagnetic shower component 
reaches its maximum is called Xmax- In hadronic interactions k/3 of the 
available energy goes into th(' electromagnetic component via Tr^'-dccays, sec 
Fig. 2. For a simple estimate only the first generation of electromagnetic 
showers is used. This approach will certainly underestimate the value of 
^max since it neglects the following subshowers. 

The first interaction occurs at an atmospheric depth Xi = In 2, where 
Xi is the interaction length of a primary proton . The latter can be 

approximated around 1 PeV by the relation 

xr--^+cigj§^ (11) 

with ^ = 68.55 g/cm^ and ( = —4.88 g/cm^. 

In the first interaction ^N^h neutral pions are produced, yielding Nch 
photons. Each photon initiates an electromagnetic cascade with the energy 
KEQ/{3Nch), developing in parallel with the others. The average multiplic- 
ity of charged particles produced in pion-nitrogen interactions [49] can be 
parameterized for energies around 1 PeV as 

with No = 55.2 and r? = 0.13. 

The depth of the shower maximum is obtained as in (1) for an electro- 
magnetic shower with an energy kEq/ (3Nch), starting after the first inter- 
action at a depth Xi, Xp,^^ = XP'"'"" ln2 + Xo\n{KEo/{3NchE^j}. Using 
(11) and (12), the expression 

XP^^^ = $ In 2 - Xo In f — + {Xq In 10 - rjXo In 10 + C In 2) lg-^° 



K PeVy ^ " ' - ■ ^ ' "pcv 

(13) 

is obtained. The elongation rate for protons is determined by the elongation 

rate for electromagnetic showers A''' = Xo In 10 and in addition by terms 
which take into account the growing multiplicity of secondary particles, as 
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Fig. 4. Average depth of the shower max- 
imum for primary photons, protons, and 
iron nuclei according to CORSIKA simu- 
lations. The lines indicate predictions ac- 
cording to (13) ( ) and the same function 

shifted up by 110 g/cm^ ( — ). 



well as the decreasing interaction length as function of energy. Taking the 
numerical parameters as described, the elongation rate A'' ~ (84.4 — 11.0 — 
3.4) g/cm^ = 70.0 g/cm^ is obtained, and one realizes that the effect of 
growing multiplicity dominates the effect of a decreasing interaction length 
by about a factor three. Evaluating also the constant term in (13) yields 
= 442.9 g/cm2 + 70.0 g/cm^ lg(i?o/PcV). 

When compared to results of full simulations, the calculated values for 
X^g^ are about 110 g/cm^ or almost 2Af~°"^ smaller than the results 
of full calculations depicted in Fig. 4 (dashed line). Presumably this is 
a consequence of neglecting the contributions of following generations of 
7r° production. However, the predicted elongation rate agrees extremely 
well with the value obtained from the CORSIKA simulations at 1 PeV 
AP = 69.9 ± 0.3 g/cm^ per decade. The solid line represents (13) shifted 
upwards by 110 g/cm^ and agrees well with the proton simulations. 

To expand the simple approach from primary protons to nuclei with 
mass number A, the superposition model is used and in (13) the energy 
is substituted by i?o M- This yields X^^^^ = X^^^~Xq In A, predicting that 
the maximum for iron induced showers should be about 150 g/ cm^ higher up 
in the atmosphere. In the full simulations, the difference is slightly smaller 
as can be inferred from the figure. 



3.1.6. Energy and Mass of the Primary Particle 

In EAS experiments the reconstructed number of electrons and muons are 
often presented in the IgTV^-lgiVg plane in order to estimate the energy 
and mass of the shower inducing particles. As an application of the simple 
Hcitler Model, lines of constant mass and energy in the N^-N^ plane are 
derived in the following to illustrate the method utilized in the experiments. 
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Fig. 5. Left: Average number of muons versus number of electrons at shower maximum 
for primary protons and iron nuclei. The data points are results of full simulations. The 
solid lines represent (14) for protons and iron nuclei, the dashed lines are equal energy 
lines according to (15). Right: Ratio of electron to muon number Ne/N/j, at shower 
maximum as function of energy for primary protons and iron nuclei. The data points are 
results of full simulations. The lines indicate (16) for the two primaries. 



To deduce lines of constant mass, (9) is transformed to obtain Eq, which, 
in turn is introduced into (6). This yields the number of muons as function 
of the number of electrons at shower maximum 

^mU i^cV^ (^^) "^^^^ ^-^^ ■ ^°-"A^e°'®** (14) 

with the exponent 5 — + 6) « 0.86. In a similar way, lines of constant 
energy are derived. A is taken from (9) and put in (6), which leads to 

AT I - — f ' p/3+^(b-i) w= ~ S 77 • 1 ^° \ ^ Ar-2-17 

E-\ a ) ^° 10 (^^p^^j A/, 

(15) 

with an exponent e = —(1 — l3)/b « —2.17. 

The constant-mass lines for protons and iron nuclei are shown in Fig. 5 
together with equal-energy lines for energies from 10^ to 10^" GeV. These 
sets of lines form a coordinate system for energy and mass in the iV^-iVg 
plane. The axis are non-perpendicular to each other. In the figure also 
results of full CORSIKA simulations for proton and iron induced showers 
are shown for fixed energies from 10^ to 3.16 • 10^° GeV in steps of half a 
decade. Taking the simplicity of the model into account the predicted lines 
agree quite well with the full simulations and they give a good illustration 
of the physics in the N^-Ne plane. 

Dividing (9) by (6) yields the electron-to-muon ratio at shower maxi- 
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mum 

It depends on the energy per nucleon Eq/A of the primary particle. This 
is the reason why the ratio N^/N/j^ is frequently used in EAS experiments 
to estimate the mass of the primary particle. If the energy is derived from 
another observable, the mass can be inferred. Predictions according to (16) 
are compared to results of full simulations for proton and iron induced 
showers in Fig. 5. The simple model predicts the calculated ratio quite well. 
The CORSIKA simulations exhibit almost a power law behavior, however, 
at high energies some flattening with respect to the predicted slope is visible. 

4. Experimental Results 

The all-particle energy spectra obtained by many experiments are com- 
piled in Fig. 1. Shown are results from direct measurements above the at- 
mosphere as well as from various air shower experiments. The individual 
measurements agree within a factor of two or three in the flux and a similar 
shape can be recognized for all experiments with a knee at energies of about 
4 PcV. Typical values for the systematic uncertainties of the absolute en- 
ergy scale for air shower experiments are about 15 to 20%. Renormalizing 
the energy scales of the individual experiments to match the all-particle 
spectrum obtained by direct measurements in the energy region up to al- 
most a PcV requires correction factors of the order of ±10% [2]. Indicating 
that the all-particle spectrum seems to be well determined. 

Due to the large fluctuations in air showers it is not possible to derive 
energy spectra for individual elements from air shower data. Therefore, 
frequently the mean mass of CRs is investigated. An often-used quantity 
to characterize the composition is the mean logarithmic mass, deflned as 
{In A) = ^^rjlnAj, r, being the relative fraction of nuclei of mass Ai. 
Investigating the ratio of the number of electrons and muons at ground 
level and the average depth of the shower maximum an increase of the mean 
logarithmic mass in the energy range around the knee could be observed 
by many experiments [1]. Such an increase is expected from consecutive 
cut-offs of the energy spectra of individual elements, starting with protons. 

A significant step forward in understanding the origin of CRs are mea- 
surements of energy spectra for individual elements or at least groups of 
elements. Up to about a PeV direct measurements have been performed 
with instruments above the atmosphere. As examples, results for primary 
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Fig. 6. Energy spectra for primary protons, helium, and iron nuclei from direct and 
indirect measurements for references see [1]. The lines indicate spectra according to the 
poly-gonato model [2]. 



protons, helium, and iron nuclei are compiled in Fig. 6. Recently, also indi- 
rect measurements of elemental groups became possible. 

With the KASCADE experiment, the problem of composition has been 
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ment, using two different models (QGSJET 01 and SIBYLL 2.1) to interpret hadronic 
interactions in the atmosphere [54] . 



approached in various ways [55]. An advanced analysis is founded on the 
measurement of the electromagnetic and muonic shower components [54] . It 
is based on the deconvolution of a two-dimensional electron muon number 
distribution. Unfolding is performed using two hadronic interaction models 
(QGSJET 01 and SIBYLL 2.1) to interpret the data. The spectra obtained 
for five elemental groups are displayed in Fig. 7. They exhibit sequential 
cut-offs in the flux for the light elements. For both models a depression 
is visible for protons around 3 to 4 PeV and at higher energies for he- 
lium nuclei. The systematic differences in flux for the spectra derived with 
QGSJET and SIBYLL amount to a factor of about two to three. The silicon 
and iron groups show a rather unexpected behavior for both models. The 
increase of the flux for both groups (QGSJET) and the early cut-off for the 
silicon group (SIBYLL) is not compatible with contemporary astrophysical 
models. The discrepancies are attributed to the fact that none of the air 
shower models is able to describe the observed data set in the whole energy 
range consistently [54]. 

The KASCADE results are compared to results of other experiments 
in Fig. 6. EAS-TOP derived spectra from the simultaneous observation of 
the electromagnetic and muonic components. HEGRA used an imaging 
Cerenkov telescope system to derive the primary proton flux [56] . The pri- 
mary proton flux has been derived from measurements of the flux of un- 
accompanied hadrons at ground level with the EAS-TOP and KASCADE 
experiments [57,58]. Spectra for protons and helium nuclei are obtained 
from emulsion chambers exposed at Mts. Fuji and Kanbala [59]. The Ti- 
bet group performs measurements with a burst detector as well as with 
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emulsion chambers and an air shower array [60,61]. 

Considering the energy range above 10 GeV, at least a qualitative pic- 
ture of the energy spectra for individual mass groups emerges: the spectra 
seem to be compatible with power laws with a cut-off at high energies. The 
spectra according to the poly gonato model [2] arc indicated in the figures 
as lines. It can be recognized that the measured values are compatible with 
cut-offs at energies proportional to the nuclear charge Ez = Z ■ 4.5 PeV. 
The lines in Fig. 1 indicate spectra for the same model. Summing up the 
flux of all elements, the all-particle flux is compatible with the flux derived 
from air shower experiments in the knee region. Above 10^ GeV the flux of 
galactic CRs is not sufficient to account for the observed all-particle spec- 
trum, and an additional, presumably extragalactic component is required. 

Energy spectra have been reconstructed with KASCADE data up to 
energies of 100 PeV. At these energies statistical errors start to dominate 
the overall error. To improve this situation, the experiment has been en- 
larged. Covering an area of 0.5 km^, 37 detector stations, containing 10 m^ 
of plastic scintillators each, have been installed to extend the original KAS- 
CADE set-up [62]. Regular measurements with this new array and the orig- 
inal KASCADE detectors, forming the KASCADE-Grande experiment, are 
performed since summer 2003 [63]. The objective is to reconstruct energy 
spectra for groups of elements up to 10^^ eV [64], covering the energy re- 
gion of the second knee, where the galactic cosmic ray spectrum is expected 
to end [65] . First analyses extend the lateral distributions of electrons and 
muons up to 600 m [66,67]. Based on one year of measurements, already 
energies close to 10^* eV are reached. It is planned to conduct an unfolding 
analysis, similar to the one described above, and reveal the energy spectra 
for groups of elements up to 10^* eV. 

A more detailed discussion of experimental results may be found else- 
where [1,26,68,69]. 

5. Conclusion and Outlook 

In the last decade the understanding of the origin of high-energy CRs has 
advanced significantly. In particular, the KASCADE experiment has shown 
that the origin of the knee in the all-particle energy spectrum is due to a 
cut-off of the light elements. A corresponding increase of the mean mass 
as function of energy in the knee region is observed by many experiments. 
Such a behavior is expected from astrophysical models, explaining the knee 
due to a finite energy reached in the acceleration process and due to leakage 
from the Galaxy. However, it has also evolved that the astrophysical inter- 
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prctation of air shower data, at present, is limited by the understanding 
of high-energy hadronic interactions in the atmosphere. Experiments hke 
KASCADE have reached the sensitivity to improve interaction models and 
corresponding analyses are under way, e.g. [70-72]. A big step forward is 
the observation of TeV-7-rays from SNRs with the expected spectral index 
7 « —2.1, thus giving an important hint to the sources of hadronic CRs. 

In the next years the KASCADE-Grande experiment and the Ice 
Cube/Ice Top experiment at the south pole [73] will measure CRs in the 
energy region of the second knee and will provide information on the mass 
composition in this region, where the galactic CR component is expected 
to end. Balloon borne experiments like ATIC, CREAM, or TRACER will 
improve the knowledge about CR propagation, by extending the energy 
spectra of individual elements to energies approaching the knee. 
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